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ABSTRACT 


The  shock-shock  interaction  problem  for  weak  in¬ 
cident  waves  impinging  on  a  supersonic  cone  has  been 
examined.  The  purpose  of  this  study  was  to  establish 
methods  for  predicting  surface  pressures  since  existing 
methods  have  been  found  inadequate  for  practical  pro¬ 
blems  with  weak  blast  waves.  It  is  shown  in  this  study 
that  one-dimensional  theory,  which  works  well  for  strong 
blast  waves,  fails  when  the  components  normal  to  the 
surface  of  the  flight  velocity  and  of  the  blast  particle 
velocity  become  comparable.  As  a  consequence,  an  axisym- 
mentric  two-dimensional  solution  method  was  developed 
using  a  "primary  wave"  approximation.  This  approach  has 
been  automated  in  the  PRIMUS  computer  code.  The  predic¬ 
tions  from  the  code  have  been  checked  against  experiment 
and  other  theory.  The  comparisons  are  presented  and  the 
agreement  is  quite  good.  The  code  was  also  used  to  pre¬ 
dict  the  cone  pressures  for  DNA  Sled  Test  experiments. 

The  results  are  presented  and  compared  with  the  experi¬ 
mental  results  and  other  theoretical  predictions. 

This  study  established  that  for  weak  blast  waves  the 
maximum  pressure  for  regular  reflection  occurs  at  the 
Mach  reflection  limit,  not  at  side-on  intercept.  The 
question  of  whether  reflection  pressures  are  even  higher 
on  the  Mach  reflection  side  of  the  boundary  was  examined, 
but  existing  prediction  methods  could  not  resolve  this 
point.  To  answer  this  question,  a  two-dimensional  finite- 
difference  code  was  also  developed  for  shock-interaction. 
The  feasibility  of  the  method  was  established,  but  suf¬ 
ficient  calculations  have  not  been  conducted  at  this 
time  to  report.  It  is  clear  that  this  finite-difference 
code  could  resolve  the  question  definitively  upon  full 
development . 
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CONCLUSIONS 


1.  The  results  from  "primary  wave"  analytical  ap¬ 
proaches  developed  for  both  one  and  two-dimensional  pro¬ 
blems  compare  well  with  exact  methods  and  with  experiment 
within  the  regular  reflection  regime. 

2.  One-dimensional  methods  have  proven  insufficient 
for  side-on  intercept  of  conical  bodies  by  weak  incident 
waves  at  supersonic  vehicle  velocities. 

3.  The  two-dimensional  analytical  method  has  been 
proven  by  comparison  with  experiment  to  give  good  esti¬ 
mates  of  reflection  pressures  on  cones  in  the  regular 
reflection  regime.  This  method  has  been  automated  as 
the  PRIMUS  (PRimary  Xnteracti°n  Method  for  Understanding 
Shocks)  computer  program. 

4.  Regular  reflection  computations  indicate  a  rapid 
rise  in  pressure  as  the  Mach  reflection  limit  is  ap¬ 
proached,  reaching  pressures  that  are  about  10%  higher 
than  in  the  side-on  situation.  These  few  computations 
suggest  that  hypersonic  similarity  parameter  correlation 
of  pressure  is  still  valid  for  side-on  intercept  in  the 
weak  blast  case  but  that  the  correlation  loses  accuracy 
near  the  regular  reflection  pressure  maximum  which  occurs 
at  the  regular-Mach  reflection  boundary. 

5.  Previously  available  analytical  and  numerical 
estimates  of  oblique  Mach  reflection  are  contradictory 
and  insufficient  to  settle  the  question:  is  there  a  maxi¬ 
mum  in  shock-shock  induced  pressure  in  the  Mach  reflection 
regime  just  before  the  transition  to  regular  reflection. 

6.  A  two-dimensional  finite-difference  method  using 
the  frozen  cells  technique  has  been  demonstrated  to  be 
feasible.  Additional  effort  and  computer  time  is  needed 
however  to  produce  useful  computations  in  both  the  Mach 
and  regular  reflection  regimes  to  settle  the  above 
question . 
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RECOMMENDATIONS 


It  is  recommended  that: 

1.  the  PRIMUS  computer  program,  that  was  developed, 
be  used  for  making  shock-shock  estimates  for  cones  and 
wedges  in  the  regular  reflection  regime, 

2.  an  effort  should  be  made  to  complete  the  two- 
dimensional  finite-difference  method  since  this  seems 

to  be  the  only  computational  method  capable  of  settling, 
once  and  for  all,  the  question  of  the  magnitude  of 
oblique  Mach  reflection  pressures,  . 

3.  should  the  completed  two-dimensional  finite- 
difference  code  demonstrate  that  the  Mach  reflection 
regime  pressure  peak  is  not  the  maximum,  it  should  be 
used  to  check  the  PRIMUS  code  in  the  regular  reflection 
regime, 

4.  if  this  check  of  the  PRIMUS  code  indicates  that 
it  needs  to  be  refined  for  a  more  exact  calculation  of 
the  complete  pressure-time  history,  an  effort  to  make 
these  refinements  should  be  undertaken  since  PRIMUS  re¬ 
quires  1/1000  (or  less)  of  the  computer  time  required  by 
the  finite-difference  method  in  the  regular  reflection 
regime, 

5.  improvement  of  these  prediction  methods  be  made 
before  or,  at  worst,  in  parallel  with  further  experi¬ 
mental  efforts  in  the  shock-shock  area. 
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INTRODUCTION 


Shock-shock  interaction  is  the  phrase  used  to  refer 
to  the  events  initiated  when  the  bow  wave  of  a  supersonic 
vehicle  encounters  an  incident  wave  from  an  external 
source,  such  as  a  blast  wave.  It  has  been  known  since 
the  early  1960's  that  such  an  interaction  may  result  in 
a  transient  pressure  rise  on  the  vehicle  just  inside  the 
incident  wave  front. 

The  concern  of  vehicle  designers  about  the  effect  of 
this  rapid  loading  and  unloading  was  originally  directed 
toward  basic  structural  integrity.  Thus,  the  mid-60's 
saw  a  number  of  development  efforts  intended  to  produce 
analytical  and/or  experimental  evaluation  of  shock-shock 
loading  for  strong  blast  waves  and  hypersonic  vehicle 
velocities.  Picatinny  Arsenal  was  actively  involved 
under  both  DAS A,  DNA,  and  AMC  sponsorship  and  produced 
analytical  studies  and  computer  codes  in  the  shock-shock 
problem  area,  (Refs  1  to  5). 

In  order  to  meet  designer's  needs,  the  SLAN  computer 
code  was  also  developed  at  Picatinny  (Ref  6) .  This  code 
is  not  an  analytical  method  but  utilizes  curve  fits  to 
the  analytical  results  developed  at  Picatinny  and  else¬ 
where.  Its  purpose  was  to  estimate  shock-shock  inter¬ 
action  loading  at  a  blunt  body  stagnation  point  and  for 
oblique  regular  and  Mach  reflection  on  cones  while  re¬ 
quiring  a  minimum  of  input  and  using  a  mimimum  of  com¬ 
puter  time. 

In  order  to  cover  the  wide  range  of  intercept  angles 
and  vehicle  geometries  cited  above,  a  great  many  assump¬ 
tions  had  to  be  made  and  hypotheses  accepted  without 
adequate  verification.  This  was  clearly  stated  and  ex¬ 
perimental  confirmation  called  for  in  the  SLAN  report. 

The  development  of  a  facility  capable  of  obtaining  the 
required  experimental  data  had  been  initiated  at  Picatinny 
Arsenal.  The  feasibility  study  had  been  completed  and  a 
pilot  facility  was  being  tested  for  performance,  (Ref  7) , 
when  all  work  on  shock-shock  interaction  was  stopped  be¬ 
cause  there  did  not  appear  to  be  any  structural  damage 
problem. 
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Shock-shock  interaction  loading  reappeared  as  a  pro¬ 
blem  in  a  different  area  of  vehicle  damage  in  1972. 

This  time,  the  concern  was  with  weak  incident  waves  and 
lower  flight  Mach  numbers.  These  conditions  invalidated 
many  of  the  assumptions  of  the  SLAN  code  -  warning  mes¬ 
sages  to  this  effect  are  printed  as  part  of  its  output. 
They  also  placed  analytically  based  numerical  approaches 
such  as  the  PASS  code  (Ref  2)  in  regimes  where  the  nu¬ 
merics  did  not  converge.  Some  outside  user's  still 
persisted  in  using  the  SLAN  code  and  in  disseminating 
the  results  without  including  the  warnings  issued  by 
the  code  itself  (for  an  example,  see  Ref  8) .  This  re¬ 
sulted  in  a  decision  to  see  what  one  and  two-dimensional 
analyses  could  be  quickly  performed  in  order  to  produce 
valid  answers  in  the  weak  wave  regime  now  of  interest. 
The  results  of  this  effort  are  presented  in  this  report. 
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DISCUSSION 


0 ne -D imensional  Methods 


The  numerical  difficulties  encountered  by  the  PASS 
code  for  weak  incident  shocks  required  that  a  different 
analytical/numerical  approach  for  incident  shock  with 
less  than  a  4  to  1  pressure  ratio  be  utilized.  At  the 
same  time,  it  appeared  desirable  to  use  as  much  already 
available  technology  as  possible. 

Picatinny  Arsenal  had  available,  a  transonic  unsteady 
finite-difference  axisymmetric  flow  field  computer  pro¬ 
gram  (Ref  9) .  Part  of  this  code  consisted  of  a  sub¬ 
routine  to  compute  which  types  of  wave  forms  will  evolve 
from  the  juxtaposition  of  two  states  of  a  gas.  Its  an¬ 
alytical  basis  is  taken  directly  from  the  analysis  begin 
ning  on  page  362  of  Landau  and  Lifshitz  (Ref  10)  except 
that  it  has  been  specialized  to  the  same  gas  on  both 
sides.  Other  subroutines  taken  from  the  transonic  pro¬ 
gram  are  used  to  interate  to  the  solution  which  matches 
the  analytical  jump  conditions  (Rankine-Hugoniot  equa¬ 
tions)  across  both  waves  and  meets  the  required  equality 
of  pressure  and  velocity  across  the  contact  surface  se¬ 
parating  them. 

The  analytical  simplifications  made  for  the  computa¬ 
tion  of  shock-shock  interaction  in  the  weak  wave  case 
were  that: 

1.  the  gas  is  ideal, 

2.  the  continuously  varying  vehicle  flow  field  prop 
erties  between  the  bow  shock  and  the  sur  face  are  re¬ 
placed  by  a  number  of  discontinuous  steps  with  the  prop 
erties  evaluated  as  the  average  of  the  end  point  values 
of  each  strip  across  the  shock  layer;  the  exceptions 
being  that  the  correct  values  are  used  at  the  shock  and 
at  the  surface  (see  Fig  1)  , 
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3 .  when  computing  the  interaction  at  each  node  only 
the  primary  wave  transmitted  toward  the  body  is  followed 
from  node  to  node,  secondary  interactions  and  waves  are 
ignored . 

The  heart  of  the  method  used  lies  in  the  second  and 
third  simplifications  above.  This  means  that  the  method 
used  is  essentially  the  method  of  Moeckel  (Ref  11)  ex¬ 
tended  to  unsteady  flow  and  it  is  also  related  to  "shock- 
expansion"  methods.  The  method  developed  will  hence¬ 
forth  be  referred  to  as  the  'primary  wave"  method.  The 
bounds  on  validity  imposed  by  the  kinds  of  simplifica¬ 
tions  made  in  the  primary  wave  method  can  only  be  deter¬ 
mined  by  comparison  with  exact  solutions  or  with  experi¬ 
ments.  The  resulting  new  one-dimensional  computer  code 
using  this  method  was  given  the  acronym  WISH  (Weak  In¬ 
teraction  of  SHocks) . 

At  the  same  time,  a  one-dimensional  "shock-tube"  un¬ 
steady  finite-difference  code  using  the  linear  method  of 
Godunov  (Ref  12)  was  modified  to  compute  the  one-dimen¬ 
sional  shock  interaction  with  the  desired  pre-interaction 
flow  field.  The  usual  difficulty  with  this  approach  is 
that  the  desired  flow  field  does  not  remain  constant  in 
time  because  it  does  not  satisfy  the  one-dimensional  flow 
equations.  It  was  found  that  this  change  occurred  more 
rapidly  than  the  shock-shock  interaction  and  wall  reflec¬ 
tion  did  for  the  cases  of  interest.  Thus,  valid  answers 
required  a  method  of  preserving  the  desired  flow  field 
for  a  sufficiently  long  time  period.  The  technique  de¬ 
veloped  here  was  a  cell-by-cell  "freezing"  of  the  desired 
flow  field  imposed  near  the  wall  end  of  the  "shock-tube." 
The  criterion  used  for  unfreezing  a  particular  cell  was  a 
1%  change  in  the  pressure  of  the  adjacent  up-stream  cell 
due  to  shock-shock  interaction. 

Both  of  the  above  methods  were  checked  by  comparing 
with  other  computations,  an  experiment,  and  with  each 
other.  The  comparison  cases  used  were: 
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1.  NOL  experiment:  flight  Mach  number  (M^  )  of  5.16, 
blast  Mach  number  (M^)  of  4.89,  head-on  intercept  of  a 
sphere,  (Ref  4) ; 

2.  PASS  code  result:  side-on  intercept  of  a  5°  cone, 

Mco=  10.0,  blast  pressure  ratio  (  )  =  21,  ideal 

air ; 


3.  Analytical:  shock  of  pressure  ratio  1.5  propag¬ 
ating  into  still  air  and  reflecting  normally  from  a 
fixed  wall. 

These  cases  were  all  computed  with  the  WISH  code  and 
with  the  one-dimensional  finite-difference  code,  which 
was  given  the  name  1-D.  A  comparison  of  results  is 
given  below: 


Table  1 


One-Dimensional  Comparison 


CASE 

NO. 

1st 

REFLECTION  PRESSURE 
(psf) 

TIME  TO  1st  WALL 
REFLECTION 
(microseconds ) 

Exp/ 

Theo 

PASS 

WISH 

1-D 

PASS 

WISH 

1-D 

1 

3,300 

(Exp) 

3,280 

3,300 

3,300 

54.4 

55.3 

57. 

2 

— 

42,360 

41,940 

41,800 

9.88 

10.0 

10. 

3 

4, 733 
(Theo) 

— 

4, 730 

4, 730 

— 

74.6 

75. 

All  results  are  seen  to  be  in  good  agreement  with 
each  other  for  both  pressure  and  time.  It  would  also 
appear  that  any  limitations  due  to  the  neglect  of 
secondary  wave  interactions  and  real  gas  effects  in 
WISH  do  not  manifest  themselves  for  nearly  equal 
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strength  shocks,  each  with  a  pressure  ratio  of  about  30 
to  1,  (Case  1)  nor  for  two  very  different  shock  strengths 
of  pressure  ratios  21  and  5  (Case  2) . 

Ihe  next  step  was  to  see  if  one-dimensional  methods 
could  produce  usable  results  for  the  side-on  intercept 
of  cones  with  weak  blast  waves  and  significant  non¬ 
parallelism  of  body  and  body  shock.  The  case  selected 
for  analysis  was  based  on  the  DNA  Shock-Shock  Test  Pro¬ 
gram  utilizing  a  rocket  sled  (Ref  8).  An  11.2°  cone  at 
Mach  5  was  intercepted  by  a  blast  wave  of  pressure  ratio 
1.6  parallel  to  the  bow  shock  of  the  cone  (16.6°  from 
the  horizontal) .  The  reflected  pressures  predicted  are 
compared  below  and  can  be  seen  to  be  in  good  agreement 
with  each  other. 


Table  2 


Sled  Test  Pressure  Comparison 


CODE 

1st  REFLECTED  PRESSURE 

WISH 

1-D 

148  psia 

143  psia 

The  questioning  of  the  validity  of  a  one-dimensional 
analysis  in  the  weak  blast  conical  case  basically  arises 
from  one  point.  Any  one -dimensional  analysis  assumes 
(among  other  things)  that  the  incident  wave,  bow  shock, 
and  body  surface  are  parallel.  The  calculations  above 
made  the  incident  shock  and  the  bow  shock  parallel  and 
let  the  cone  surface  become  non-parallel.  Next  the  WISH 
code  was  used  to  compute  the  case  where  the  incident 
shock  and  the  cone  surface  were  parallel  but  the  bow 
shock  was  not.  The  result  was  a  wall  reflection  pressure 
of  95  pounds  per  square  inch.  This  large  difference  in¬ 
dicates  that  one-dimensional  analysis  is  not  adequate. 
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The  question  then  becomes  why  was  a  one-dimensional 
analysis  adequate  for  stronger  incident  waves  (pressure 
ratios  of  5  to  201)  and  thin  vehicle  shock  layers  (Ref 4) ? 
The  answer  seems  to  lie  in  the  relative  values  of  the 
blast  particle  velocity  and  the  component  of  cone  flight 
velocity  normal  to  the  blast  front  anci,  hence,  in  momen¬ 
tum  and  energy.  When,  for  example.  Check  Case  #2  and 
the  Sled  Test  input  conditions  to  the  WISH  code  are  com¬ 
pared,  the  numbers  are  those  given  below. 

Table  3 

SIDE-QN  INPUT  CONDITIONS  COMPARISONS 


CASE 

CHECK  CASE  #2 

DNA  SLED  TEST 

Vehicle 

Mach  No . 

10.0 

5.0 

Blast 

Pressure  Ratio 

21.0 

1 . 6 

Blast  Particle 
Velocity 

3274  fps 

388  fps 

Cone 

Ha If -Angle 

5.0° 

11.2° 

Cone  Shock 

Angle 

7.7° 

16.6° 

Flight  Velocity 

J_  Bow  Shock 

1296  fps 

1588  fps 

Flight  Velocity 

JL  Cone  Surface 

843  fps 

1080  fps 

It  is  clear  that  the  major  difference  between  the 
cases  is  that  in  the  strong  wave  case  the  blast  particle 
velocity  is  70  to  80%  of  the  sum  of  normal  vehicle 
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velocity  and  blast  particle  velocity,  while  in  the  weak 
wave  case  the  proportions  are  just  about  reversed  and 
the  vehicle  velocity  component  dominates.  In  the  latter 
situation  the  choice  of  whether  the  blast  wave  is  paral¬ 
lel  to  the  body  surface  or  parallel  to  the  body  shock 
has  a  substantial  effect  upon  the  resulting  total  veloc¬ 
ity  in  body-fixed  coordinates.  Conversely,  in  the  strong 
incident  shock  case,  the  total  velocity  is  relatively 
insensitive  to  whether  the  blast  is  aligned  with  the  body 
surface  or  with  the  body  shock  since  so  much  of  the  total 
velocity  is  due  to  the  blast  particle  velocity.  It  is 
also  clear  that  the  thicker  the  shock  layer,  the  greater 
the  difference  between  the  component  of  the  flight  veloc¬ 
ity  normal  to  the  body  surface  and  the  component  normal 
to  the  body  shock. 

As  a  result  the  use  of  one-dimensional  analyses  for 
the  weak  wave  case,  which  is  dominated  by  a  two-dimen¬ 
sional  flow,  is  not  appropriate. 

Two-Dimensional  Methods  (Planar  and  Axisymmetric) 

It  is  well  established  from  experiment  (Ref  13  to  15) 
that  the  initial  interaction  between  a  conical  shock  and 
an  incident  wave  can  be  treated  as  a  planar  two-dimen¬ 
sional  problem.  This  interaction  can  be  solved  exactly 
by  the  SWIVEL  code  available  at  Picatinny  Arsenal.  The 
SWIVEL  code  is  also  capable  of  computing  the  regular  re¬ 
flection  pressure  and  duration  of  the  transient  due  to 
the  transmitted  incident  wave  when  the  wave  propagates 
through  a  constant  property  layer  after  the  shock-shock 
interaction,  i.e.,  the  case  of  a  wedge,  (See  Fig  2  for 
the  basic  wave  diagram  and  designation  of  regions) . 

The  non-uniform  shock  layer  of  a  cone  was  treated 
by  the  same  basic  method  of  using  a  step-wise  approxima¬ 
tion  of  a  non-uniform  shock  layer  with  which  the  trans¬ 
mitted  incident  wave  interacts  and  ignoring  the  secondary 
interactions.  The  SWIVEL  code's  analytical  basis  is  the 
computation  of  the  4-shock  intersection  that  can  result 
when  2  co-planar  waves  interact.  Thus  in  the  new  analy¬ 
sis,  the  isen tropic  compression  between  shock  and  surface 
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on  a  cone  was  modeled  by  a  series  of  weak  shock  compres¬ 
sions  since,  as  is  well  known,  weak  shocks  are  very 
nearly  isentropic. 

During  the  development  of  the  new  computer  code,  it 
was  found  that  a  single  interior  shock  compression  to 
the  cone  surface  conditions  gave  essentially  the  same 
reflection  pressure  as  a  multi-shock  compression.  It 
was  also  found  that  for  the  cone  fL  ight  conditions  of 
interest,  that  the  entropy  jump  across  a  single  interior 
shock  was  still  small  enough  that  the  flow  was  nearly 
isentropic.  Therefore,  in  the  interest  of  simplicity 
of  coding,  getting  results  sooner,  and  ease  of  setting 
up  cases  a  single  shock  compression  was  used.  The  new 
computer  code  was  given  the  acronym  PRIMUS  (PRimary  Inter¬ 
action  Method  of  Understanding  Shocks) .  PRIMUS  also 
retains  the  wedge  case  capability  of  the  SWIVEL  code. 

Another  change  that  had  to  be  made  in  SWIVEL  to  pro¬ 
duce  PRIMUS  was  the  addition  of  the  ability  to  compute 
the  quasi-steady  conditions  on  a  cone  at  an  angle  of 
attack.  The  wedge  conditions  can  be,  and  still  are,  com¬ 
puted  in  closed  form  but  the  cone  conditions  at  angle  of 
attack  can  not.  Windward  ray  cone  surface  pressures  at 
angle  of  attack  are  estimated  from  experimental  data  col¬ 
lected  and  correlated  by  one  of  the  authors  in  the  form 
of  the  ratio  of  actual  pressure  coefficient  to  the  New¬ 
tonian  pressure  coefficient  vs  the  hypersonic  similarity 
parameter  defined  by  M^  sin  (0c  +  oO»  where  M^  is  the 
flight  Mach  number,  0C  is  the  cone  semi-angle,  and  o<  is 
the  angle  of  attack.  Other  properties  required  were 
taken  from  Sims'  tabulation  for  cones  at  small  angles  of 
attack  (Ref  16)  and  curve  fitted.  The  curve  fits  are 
restricted  to  the  range  of  parameters  currently  expected 
to  be  required  and  they  have  not  been  checked  for  ac¬ 
curacy  at  flight  Mach  numbers  less  than  3,  or  cone  half¬ 
angles  greater  than  15°.  Furthermore,  Sims'  angle  of 
attack  data  is  only  valid  at  small  angles  of  attack. 

The  original  SWIVEL  code  estimated  the  duration  of 
the  pressure  transient  by  assuming  that:  the  reflected 


12 


wave  was  not  deflected  or  altered  in  strength  by  its  in¬ 
teraction  with  the  contact  surface,  the  expansion  wave 
resulting  from  the  interaction  of  the  reflected  wave  and 
the  deflected  blow  wave  also  propaged  back  to  the  body 
through  uniform  conditions,  and  that  this  expansion  was 
the  primary  pressure  relief  mechanism.  The  duration  was 
defined  so  that  the  decay  time  due  to  the  width  of  the 
expansion  fan  was  excluded.  This  analysis  has  not  been 
altered  at  this  time  in  PRIMUS.  Thus,  the  times  pre¬ 
dicted  for  the  conical  case  have  been  computed  using  the 
same  basic  assumptions.  Of  course  the  uniform  conditions 
used  are  the  true  cone  surface  conditions.  A  desirable 
future  effort  would  be  to  incorporate  the  complete  non- 
uniform  layer  including  interaction  with  the  contact  sur¬ 
face  . 

In  parallel  with  this  development,  an  effort  was  made 
to  modify  a  two-dimensional,  unsteady  finite-difference 
blunt  body  flow  field  program  (Ref  17) ,  developed  with 
DNA  funding,  to  compute  shock-shock  interaction  of  a 
plane  shock  wave  with  a  conical  flow  field.  This  ap¬ 
proach  allows  the  realistic  approximation  of  two-dimen¬ 
sional  effects  by  allowing  for  non-parallel  body  surface, 
bow  shock  and  blast  wave  and  for  the  non-uniform  shock 
layer,  which  is  the  main  flow  difference  between  wedges 
and  cones.  This  is  modeled  by  the  "freezing"  of  the  im¬ 
posed  conical  flow  just  as  has  been  discussed  for  the  one¬ 
dimensional  finite  difference  code. 

This  approach  is  the  only  one  to  date  which  offers 
hope  for  answering  certain  questions  that  have  been  asked 
in  the  technical  community  as  to  where  the  maximum  wall 
pressure  occurs  as  a  function  of  intercept  angle:  is  there 
a  maximum  on  the  Mach  reflection  side  of  the  transition 
intercept  angle  from  Mach  to  regular  reflection  and  is 
it  greater  than  the  rise  in  pressure  on  the  regular  re¬ 
flection  side  of  the  transition.  The  finite  difference 
method  is  expected  to  compute  both  with  equal  facility. 

Its  one  serious  drawback  would  seem  to  be  long  running 
times  (about  one  hour  on  a  CDC  6500  which  would  make  it 
a  candidate  for  answering  certain  questions  about  pheno- 
meonology  rather  than  a  code  to  be  used  for  each  case  of 
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interest.  The  "frozen  flow  field"  used  in  the  one¬ 
dimensional  finite-difference  code  and  now  incorporated 
into  the  two-dimensional  code  does  not  affect  the  pro¬ 
pagation  rate  of  discontinuities  nor  reflected  pressures. 
All  the  unsteady  finite-difference  methods  of  lower  di¬ 
mensionality  than  the  real  problem  do  share  one  dif¬ 
ficulty:  they  compute  the  final  quasi-steady-state 
solution  of  the  lower -dimensionality  problem  rather  than 
that  of  the  problem  they  are  being  used  to  approximate. 
This  is  not  an  insuperable  problem  and  it  can  be  resolved 
by  the  same  non-numer ical  method  used  in  the  PASS  code; 
one  estimates  the  width  of  the  final  expansion  wave  and 
its  strength  by  simple  closed  form  methods  (page  78  of 
Ref  18)  and  patches  this  onto  the  finite  difference  solu¬ 
tion  beginning  at  the  first  predicted  pressure  relief 
at  the  body  surface. 

This  study  was  carried  to  the  point  of  showing  fea¬ 
sibility.  The  code  was  modified  to  planar  two-dimen¬ 
sional  form,  the  boundary  conditions  suitably  modified 
and  initial  conditions  for  the  shock-shock  interaction 
problem  coded  and  debugged.  A  true  wedge  case  for  re¬ 
gular  reflection  was  set  up  (as  the  correct  result  was 
known)  and  the  computation  begun.  The  initial  inter¬ 
action  of  the  bow  and  blast  waves  proceeded  normally  and 
the  transmitted  blast  wave  was  propagated  toward  the 
body  surface;  then  the  body  surface  pressure  near  the 
nose  began  to  rise  as  it  should.  Unfortunately,  at  this 
time  development  work  had  to  cease.  The  method  appears 
to  be  verified  as  to  feasibility  and  a  small  effort 
would  determine  its  accuracy  for  wedge  regular  reflection, 
cone  side-on  at  high  Mach  number,  and  Mach  reflection 
for  cones  or  wedges,  both  head-on  and  oblique. 

Mach  Reflection 

The  question  has  arisen,  provoked  by  some  predictions 
of  Gardner  (Ref  20)  and  two  experimental  points  in  Ref  19, 
as  to  where  the  maximum  shock-shock  interaction  surface 
pressure  as  a  function  of  intercept  angle  actually  occurs. 
Gardner  at  McDonnell-Douglas  suggests  that  it  lies  at  the 
boundary  between  Mach  and  regular  reflection  on  the  Mach 
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reflection  side  and  may  exceed  any  regular  reflection 
pressure.  To  illustrate  this,  the  results  from  Ref  20 
has  been  redrawn  as  Fig  3  of  this  report. 

The  experimental  pressure  ratios  in  Ref  19  are  for 
a  9°  half-angle  cone  at  Mach  3.1  and  a  blast  pressure 
ratio  of  5.  The  PRIMUS  code  predicts  the  result  for 
HARTS  Round  7  (40°  intercept  angle)  to  within  the  in¬ 

dicated  uncertainty  and  .predicts  the  limit  for  regular 
reflection  to  be  between  39.5°  and  39.0°.  It  is  HARTS 
Round  9,  which  was  supposedly  at  a  41°  angle  of  inter¬ 
cept,  which  gave  an  80%  higher  value  of  reflected  pres¬ 
sure.  Since  the  angle  of  attack  of  the  cone  and  angle 
of  intercept  of  the  blast  are  both  subject  to  some  un¬ 
certainty,  this  may  be  evidence  for  a  maximum  in  pres¬ 
sure  just  on  the  Mach  ref lection  side  of  the  limit  if 
the  combined  effect  was  to  reduce  the  effective  angle 
of  intercept  below  about  39°.  On  the  other  hand,  it 
may  be  evidence  for  an  error  in  data  acquisition  or  re¬ 
duction. 

A  brief  survey  of  existing  analytical  work  that 
applies  to  oblique  Mach  reflection  reveals  very  little 
of  direct  applicability  to  the  cases  of  interest. 
Blankenship  and  Busemann  (Ref  21)  present  some  calcula¬ 
tions  for  cones  with  very  small  intercept  angles  (1°  or 
less  away  from  head-on) .  These  results  apply  to  a  wave 
pattern  which  occurs  only  at  very  high  vehicle  Mach 
numbers  (M<?0»  10)  for  the  blast  strengths  of  interest. 
These  calculations  can  only  give  an  initial  slope,  which 
certainly  does  not  remain  constant,  for  the  variation 
with  intercept  angle  while  the  intercept  angles  of  inter¬ 
est  are  20  to  30  degrees.  Most  fundamentally,  the  as¬ 
sumed  wave  pattern  does  not  apply  to  the  flight  and  blast 
Mach  numbers  of  current  interest.. 

The  wedge  analysis  and  results  by  Inger  (Ref  22) 
appears,  at  first  glance,  to  be  applicable.  His  assump¬ 
tion  of  a  weak  blast  wave  is  consistent  with  the  current 
problem  and  the  same  method  is  used  for  intercept  angles 
from  0°  (head-on)  to  60°.  The  assumption  of  a  wedge  need 


15 


not  disqualify  this  method  for  useful  cone  computations, 
(see  the  comparison  in  Fig  4).  However,  Inger's  method 
also  requires  a  strong  bow  wave  and  he  appears  unwill¬ 
ing  to  extend  his  method  below  a  value  of  the  hypersonic 
similarity  parameter  of  about  2  (see  Fig  7  of  Ref  22). 
Upon  closer  examination,  one  finds  that  the  restrictions 
are  even  more  severe.  Not  only  must  the  hypersonic 
similarity  parameter  be  large  but  this  must  be  achieved 
while  keeping  the  wedge  and  wedge  shock  essentially 
parallel.  This  implies  small  wedge  angles  and  high  Mach 
numbers.  Thus  the  range  of  proposed  validity  seems  to 
be  at  too  high  Mach  numbers  and  for  too  slender  wedges 
to  be  of  use  in  the  current  problem  area.  Nevertheless, 
one  might  hope  to  at  least  answer  the  important  question 
of  the  trend  of  peak  pressure  with  intercept  angle  as 
it  increases  from  head-on.  Inger's  method  shows  a 
smooth  steady  increase  in  peak  pressure  with  increasing 
obliquity  of  intercept. 

The  stated  range  of  validity  is  that  the  original 
wedge  shock  angle  (  @  )  minus  the  wedge  angle  (  S  ) , 

(all  in  radians,  see  also  Fig  2)  is  much  less  than  1. 
However,  the  results  seem  erroneous,  at  least  for  large 
intercept  angles,  even  in  this  proposed  range  of  vali¬ 
dity.  The  PRIMUS  code  computes  the  exact  solutions  for 
shock-shock  interactions  on  wedges  which  result  in  reg¬ 
ular  reflection  from  the  surface.  A  check  case  of  a  6° 
wedge  at  Mach  20  (for  which  the  angular  difference  is 
0.04  radians)  with  a  blast  wave  of  Mach  number  1.2  was 
computed  with  PRIMUS  for  intercept  angles  from  60°  to 
22°.  Results  for  this  case  were  also  computed  using 
the  correlations  plots  in  Ref  22.  A  comparison  of  re¬ 
flection  pressure/ambient  pressure  is  shown  below  and 
plotted  in  Fig  5. 
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Table  4 


Oblique  Intercept  Comparison  for  Wedges 


INTERCEPT  ANGLE 

prefl) 

prefl) 

(0°  =  head-on) 

\  P  co  /  Inger 

|pw  /  PRIMUS 
(Exact) 

f— ' 

o 

0 

12.60 

Mach  Reflection 

20° 

13.10 

Mach  Reflection 

25° 

13 .38 

16.11 

30° 

13.68 

15.89 

0 

o 

14.41 

16.01 

50° 

15.45 

16.24 

60° 

17.11 

16.45 

Since  the  same  pre-intercept  and  quasi-steady  pres¬ 
sures  were  used  for  both  methods  for  each  case,  it  seems 
that  the  difference  must  be  attributed  to  the  value  of 
maximum  pressure  computed  by  Inger's  method.  That 
method  does  not  predict  the  minimum  in  pressure  followed 
by  a  strong  rise  with  decreasing  intercept  angle  just 
before  the  Mach  reflection  limit  that  is  characteristic 
of  weak  incident  wave  regular  reflection.  Either  Inger's 
method  cannot  be  applied  to  as  large  intercept  angles  as 
he  believed  in  1966  or  0.04  is  not  «  1  for  his  method 
at  =  20.  Since  the  small  pertubation  theory  used 

is  only  valid  in  the  limit  as  6— ►  0,  this  could  be  the 
case.  Whether  or  not  it  is  adequate  for  small  inter¬ 
cept  angles  (Mach  reflection)  is  not  well  established  at 
this  time. 
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A  rather  risky  extrapolation  of  cone  experimental 
results  in  Ref  23  to  lower  blast  strengths  and  assuming 
that  the  hypersonic  similarity  parameter  is  a  correlat¬ 
ing  quantity  indicates  that  Inger 1 s  results  for  head-on 
intercept  in  this  case  are  in  good  agreement  (see  Fig  5). 
If  this  is  a  valid  check,  and  if  Inger  predicts  the  cor¬ 
rect  trend,  it  indicates  that  there  is  no  reflection 
pressure  maximum  due  to  Mach  reflection.  Any  rise  in 
pressure  near  the  Mach-regular  reflection  limit  would 
then  be  due  to  the  rise  which  occurs  in  the  regular 
regime.  However,  there  are  too  many  unconfirmed  as¬ 
sumptions  to  consider  the  matter  settled. 

The  work  of  Smyrl  (Ref  24)  is  similar  to  that  of 
Inger  and  is  applicable  to  lower  flight  Mach  numbers 
but  it  is  a  linearized  analysis  which  requires  that  the 
wedge  angle  (or  wedge  angle  plus  the  intercept  angle) 
be  small  enough  that  linearization  is  valid  at  the  de¬ 
sired  Mach  number.  Thus,  this  method  can  only  give  the 
initial  slope  of  the  pressure-intercept  angle  function 
for  wedges  that  are  intercepted  obliquely. 

RESULTS 


The  one -dimensional  "primary  wave"  (WISH  code)  and 
finite-difference  methods  developed  (1-D  code)  were 
found  to  be  in  good  agreement  with  each  other  and  with 
a  number  of  other  exact  and  experimental  results  for  one¬ 
dimensional  flow.  One -dimensional  computations  of  the 
DNA  Sled  Test  conditions  at  a  Mach  number  of  5,  an  11.2° 
half-angle  cone  with  a  16.6°  bow  wave  angle,  and  a  blast 
pressure  ratio  of  1.6  indicated  that  there  was  a  55%  in¬ 
crease  in  reflected  surface  pressure  in  going  from  a 
blast  wave  parallel  to  the  cone  surface  to  a  blast  wave 
parallel  to  the  bow  shock.  Reasons  for  this  difference 
are  detailed  in  the  Discussion  and  the  essential  result 
is  that  the  one -dimensional  approximation  breaks  down 
at  some  difference  in  angle  when  the  three  surfaces 
(blast  wave,  bow  wave,  and  body)  are  not  parallel.  The 
breakdown  occurs  at  small  angles  when  the  vehicle  velocity 
norma 1  to  the  blast  wave  dominates  the  blast  particle 
velocity  in  the  sum  of  the  two  terms. 
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The  two-dimensional  "primary  wave"  analytical  ap¬ 
proximation  was  automated  in  a  form  appropriate  for 
cones  in  the  PRIMUS  computer  program.  Suitable  experi¬ 
mental  or  independent  analytical  methods  to  check  against 
are  harder  to  obtain  for th e oblique  intercept  of  cones. 

The  most  suitable  predictions  available  were  those  by 
Gardner  and  Lenertz  (MDAC  and  MMC,  respectively,  on 
Fig  2-15  of  Ref  8)  for  the  DNA  Sled  Test.  The  PRIMUS 
results  are  in  good  agreement  with  Gardner  for  pressure 
and  in  fair  agreement  for  duration  in  this  case  with  the 
blast  front  nearly  parallel  to  the  cone  surface.  This 
is  shown  in  Figure  6.  Lenertz  does  not  predict  durations 
and  the  pressure  predicted,  is  lower  than  the  other  methods. 
The  peak  pressure  predicted  by  PRIMUS  is  73  to  79  psia, 
depending  on  the  experimental  uncertainty  in  blast 
strength  and  orientation. 

The  experimental  results  for  the  Sled  Test  are  re¬ 
ported  in  Reference  8.  The  data  in  which  the  most  con¬ 
fidence  is  placed  is  from  the  test  designated  4B-B3 . 
However,  the  fact  that  the  high  temperature  pressure 
gauges  used  had  too  slow  a  response  time  and  the  regular 
gauges  were  subjected  to  a  high  temperature  environment 
had  necessitated  a  great  deal  of  processing  of  the  raw 
data.  The  basic  technique  was  to  calibrate  the  gauges 
with  a  shock-tube  and  use  convolution  techniques  to  pre¬ 
dict  the  actual  pressures  sensed  during  the  rise  time 
of  the  slower  response  gauges.  Figure  2-15  of  Ref  8 
shows  both  the  preliminary  estimate  band  and  the  final 
estimate  of  the  shock-shock  reflection  pressure  after 
processing  the  raw  data.  The  final  estimate  results  in 
a  peak  reflected  pressure  of  about  100  psia.  This  is 
3  2%  above  the  me  an  estimates  of  any  of  the  two-dimensional 
theories  including  PRIMUS.  The  final  reduced  data  from 
the  experiment,  as  taken  from  Ref  8,  are  compared  with 
the  theoretical  predictions  on  Fig  6.  It  should  be 
noted  however  that  the  gauge  in  question  read  the  theo- 
etical  pre-interaction  surface  pressure  (33  psia)  ex¬ 
actly  while  it  read  23%  high  for  the  quasi-steady  post¬ 
interaction  pressure.  This  suggests  that  the  gauge  had 
changed  calibration  during  the  shock-shock  transient. 

The  PRIMUS  predictions  and  an  analytical  approximation 
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from  Reference  14  are  compared  next  with  a  different 
experiment  (Ref  14)  which  did  not  depend  upon  pressure 
gauges.  Since  good  agreement  was  achieved  with  this 
other  experiment,  and  several  theories  are  in  essential 
agreement  for  the  DNA  Sled  Test,  it  would  seem  that  the 
gauge  response  problem  is  the  major  cause  of  the  dis¬ 
agreement  between  theory  and  experiment  in  that  test. 
Certainly,  the  experiment  was  not  as  definitive  as  one 
would  have  liked.  It  is  understood  that  a  Phase  II  ex¬ 
periment  will  be  conducted  in  which  it  is  hoped  to  eli¬ 
minate  some  of  the  problems  mentioned. 

The  only  oblique  regular  reflection  on  cones  experi¬ 
mental  results  are  those  of  Ref  14  where  pressures  are 
predicted  from  photographs  of  shock  patterns  observed 
on  cone  models  launched  through  an  oblique  shock  at  40° 
from  the  normal  to  the  cone  axis.  The  incident  shock 
strengths  tested  are  higher  than  those  of  interest  to 
the  current  problem  but  the  cone  angles  and  cone  Mach 
numbers  are  in  the  desired  range. 

The  PRIMUS  predictions  and  these  experimental  re¬ 
sults  for  the  ratio  of  reflected  pressure  to  initial  cone 
pressure  are  compared  in  Figures  7  to  10.  Only  the  mean 
value  of  the  experimental  shock  wave  Mach  numbers  actu¬ 
ally  achieved  was  used  at  each  nominal  experimental  value 
for  the  incident  wave  strength  in  the  PRIMUS  computations. 
Quite  good  agreement  was  found  for  the  variation  of  the 
ratio  of  reflected  to  initial  pressure  with  shock  wave 
Mach  number  for  both  9°  and  15°  half-angle  cones  (see 
Fig  7  and  8) .  The  PRIMUS  predictions  of  the  variation 
with  cone  Mach  number  shows  a  somewhat  stronger  depend¬ 
ence  on  cone  Mach  number  than  the  experimental  data  in¬ 
dicates  for  both  9°  and  15°  cones,  at  the  higher  blast 
strengths.  As  the  blast  strength  is  reduced,  the  agree¬ 
ment  with  experiment  becomes  better  and  better  and  is  ex¬ 
cellent  at  the  lowest  strength  tested  and  should  be  ex¬ 
cellent  at  the  even  lower  strengths  of  interest,  (see 
also  Fig  9  and  10) * 
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The  Two-dimensional  theory  used  in  Ref  14  computes 
the  2 -shock  interaction  exactly  but  ignores  the  effect 
of  the  non-uniform  shock  layer  on  the  strength  and  orient¬ 
ation  of  the  transmitted  shock  and  then  computes  the  re¬ 
flection  at  the  surface  using  actual  cone  surface  condi¬ 
tions.  The  primary  wave  approach  in  PRIMUS  does  con¬ 
sider  the  non-uniformity  of  the  shock  layer.  The  im¬ 
provement  made  by  accounting  for  the  change  in  trans¬ 
mitted  wave  strength  and  inclination  is  discussed  below. 
The  predictions  of  the  theory  of' Ref  14  are  always  higher 
than  the  PRIMUS  code  predictions.  The  PRIMUS  predictions 
are  clearly  in  better  agreement  with  experiment  for  the 
9°  cones  (Fig  9);  while  in  the  case  of  the  15°  cones, 
the  data  scatter  is  too  great  to  permit  a  choice  between 
the  two  theories  (Fig  10) .  The  pressure  trends  with 
shock  wave  Mach  number  of  the  two  methods  can  be  seen  by 
comparing  Fig  7  with  11  and  Fig  8  with  Fig  12.  Again  the 
method  used  in  PRIMUS  compares  more  favorably  with  ex¬ 
periment  for  the  9°  cones  and  also  appears  superior  for 
the  15°  cones. 

The  PRIMUS  code  achieved  this  agreement  even  though 
various  messages -warning  of  the  violation  of  assumptions- 
were  printed  out  for  some  of  the  required  cases.  How¬ 
ever,  these  were  related  to  the  small  angles  of  attack 
limit  inherent  in  using  Sims'  results  (Ref  16)  for  the 
quasi-steady  state  or  to  the  duration  calculation,  which 
is  known  to  be  an  approximation.  It  is  hoped  that  this 
will  not  lead  to  the  wholesale  neglect  of  such  messages 
by  users  as  it  did  with  the  SLAN  code. 

The  two-dimensional  finite  difference  computer  code 
using  a  method  of  "freezing"  the  imposed  preinteraction 
flow  field  was  brought  to  the  point  of  proving  feasibil¬ 
ity  of  the  basic  method  in  two  dimensions.  The  necessity 
of  stopping  work  prevented  ch  ecking  it  against  exact  solu¬ 
tions  for  wedges  and  comparing  it  with  the  PRIMUS  results 
for  cones.  This  computer  code  is  the  only  one  which  cur¬ 
rently  has  the  potential  of  computing  the  oblique  Mach 
reflection  case  and,  hence,  resolving  the  uncertainties 
in  this  regime  dealt  with  in  detail  in  the  Discussion . 
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In  order  to  permit  the  ready  computation  of  oblique 
intercept  of  cones  and  wedges  in  the  regular  reflection 
regime.  Appendix  A  of  this  report  contains  a  brief 
user's  guide  to  input  and  output  of  the  PRIMUS  code,  a 
listing  of  PRIMUS,  and  a  wedge  and  a  cone  check  case. 

The  PRIMUS  code  has  been  used  to  compute  shock- 
shock  results  for  a  cone  of  semi-angle  11.2°  and  a 
blast  Mach  number  of  1.234  at  intercept  angles  from 
90°  down  to  the  Mach  reflection  boundary  for  fliqht 
Mach  numbers  of  3,  5,  and  7.  The  variation  of  the  maxi¬ 
mum  in  surface  pressure  over  ambient  pressure  with  in¬ 
tercept  angle  and  flight  Mach  number  is  shown  in  Fig  13 . 
The  maximum  reflected  pressure  ratio  at  each  Mach  number 
is  plotted  versus  flight  Mach  number  on  Fig  14.  The 
pressure  ratio  at  an  intercept  angle  of  78.8°  (blast 
wave  parallel  to  cone  surface)  is  also  plotted  against 
Mach  number  on  the  same  figure.  This  shows  that  the 
maximum  pressure  in  the  regular  reflection  regime  always 
occurred  at  the  boundary  between  Mach  and  regular  re¬ 
flection,  not  at  the  "side-on"  conditions,  for  weak 
blast  waves.  At  the  same  blast  strength  and  cone  semi¬ 
angle,  the  regular  reflection  boundary  moves  to  lower 
intercept  angles  (more  nearly  head-on)  as  the  flight 
Mach  number  increases.  This  is  due  to  the  stronger  bow 
wave  resulting  in  greater  weakening  of  the  transmitted 
wave.  At  the  same  time,  however,  the  reflected  pressure 
over  ambient  pressure  ratio  rises  with  Mach  number  since 
the  surface  pressure  rises  more  strongly  with  flight 
Mach  number . 

The  rapid  rise  in  reflection  pressure  with  small 
changes  in  intercept  angle  near  the  Mach-regular  reflec¬ 
tion  boundary  can  result  in  significant  errors  if  the 
computational  angular  increment  is  too  large  in  this 
region.  For  example  in  the  above  computations  at  M^  =  5, 
if  calculations  are  made  at  one  degree  intervals,  a 
ratio  of  reflection  pressure  to  ambient  of  6.02  at 

=  25°  and  Mach  reflection  at  24°  is  obtained.  But 
when  increments  of  0.1°  are  used  the  pressure  ratio  be¬ 
comes  6.36  at  24.2°,  with  Mach  reflection  at  24.1°. 


22 


When  increments  of  0.01°  are  used,  the  maximum  pressure 
ratio  increases  still  further  to  6.44  at  24.18°  with 
Mach  reflection  at  24.17°  (see  also  Fig  15). 

Earlier  side-on  results  for  strong  blast  waves  had 
correlated  well  with  the  hypersonic  similarity  para¬ 
meter,  (Ref  4).  In  order  to  see  if  this  remained  true 
for  weak  blasts,  the  Sled  Test  case  and  a  case  at  twice 
the  Mach  number  but  the  same  value  of  the  hypersonic 
similarity  parameter  were  computed  with  PRIMUS.  The  re¬ 
sults  are  compared  on  Fig  15  for  the  same  blast  strength 
at  various  intercept  angles.  The  side-on  ratios  of  re¬ 
flection  pressure  to  ambient  are  nearly  equal,  however, 
the  Mach  10  results  first  undershoot  then  overshoot  the 
Mach  5  results  as  the  intercept  angle  decreases.  The 
maximum  error,  at  the  same  anale,  occurs  near  the  Mach 
reflection  limits  and  is  about  10%. 
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APPENDIX  A 


PRIMUS  INPUT 

DESCRIPTION 

COLUMN 

NAME 

FORMAT 

EXPLANATION 

1-10 

Ml 

F 

Flight  Mach  number 

11-20 

XIB 

F 

Blast  pressure/ambient 
pressure^ 

21-30 

BETA 

F 

Blast  intercept  angle 
^qO  =  head-on) 

31-40 

DELTA 

F 

Wedge  case  =  wedge  angle; 
deg.  Need  not  be  specified 
if  XI 1  is  specified.  Cone 
Case  =  flow  deflection  angle 
at  conical  shock,  deg.  need 
not  be  specified  if  XIS  is 
specified . 

41-50 

XI 1 

F 

Pressure  ratio  across  wedge 
shock.  Need  not  be  specified 
if  DELTA  is  specified  for 
wedge.  Surface  pressure 
ratio  for  cone  case. 

51-60 

CONE 

F 

Cone  half-angle,  deg.  Blank 
for  wedge  case 

61-70 

XIS 

F 

Cone  shock  pressure  ratio. 
Need  not  be  specified  if 
DELTA,  conical  flow  deflec¬ 
tion  angle,  is  specified. 
Blank  for  wedge  case. 

^•Analysis  currently  scaled  to  ambient  pressure  =  1,  am¬ 
bient  sound  speed  =  1,  ambient  density  =  ratio  of  specific 
hearts  =  1.4. 
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PRIMUS  KEY  OUTPUT  QUANTITIES1 

MSUBl 

flight  Mach  number 

XISUBB 

blast  pressure/ambient  pressure 

BETA 

incident  wave  intercept  angle  (0°  =  head-on) 

deg. 

GAMMA 

ratio  of  specific  heats  (constant) 

DELTA 

wedge  angle,  or  flow  deflection  angle  at 
conical  shock,  deg. 

ASUBl 

ambient  sound  speed  (scaled  to  1) 

PSUBl 

ambient  pressure  (scaled  to  1) 

MSUBB 

blast  wave  Mach  number 

LAMBDA 

final  state  angle  of  attack,  deg. 

THETAF 

final  state  wedge  shock  angle,  deg. 

XISUBWF 

final  state  wedge  pressure/blast  pressure 

PSUB3F 

final  state  wedge  pressure/ambient  pressure 

P6/P3F 

final  state  wedge  reflection  pressure/ambient 
pressure 

PSUBF 

cone  or  wedge  reflection  pressure/ambient 
pressure 

MSUBF 

final  state  cone  surface  Mach  number 

THETAC 

final  state  cone  shock  angle,  windward  ray,  deg 

ALPHAC 

final  state  cone  angle  of  attack,  deg. 

CONANGL 

cone  half  angle,  deg. 

28 


pc/pl 

cone  initial  surface  pressure/ambient  pressure 

P6/PC 

cone  reflection  pressure/initial  cone  surface 
pressure 

TAU 

cone  or  wedge  dimensionless  duration  of  re¬ 
flection  pressure,  ta^/b,  where  t  is  the  actual 

time  and  b  is  the  axial  distance  from  the  wedge 

or  cone  tip  to  the  point  of  interest.^ 


i)See  also  Fig  2  for  designation  of  regions  and  points. 
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LISTING  OF  PRIMUS 
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V 


CHECK  CASES  FOR  PRIMUS 


1.  Wedge 

Flight  Mach  Number  =  10.,  Wedge  Angle=ll . 536° 

Blast  Pressure  Ratio  =  1.513,  Intercept  Angle=70° 


INPUT  CARD 


Flight  Mach  Number  =5.,  Cone  Semi-Angle=11^2° 

Blast  Pressure  Ratio  =  1.64,  Intercept  Angle=70 


INPUT  CARD 
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NON-DIMENSIONAL  DENSITY 


DENSITY  VARIATION  ACROSS  SHOCK  LAYER 
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FIG. I 


BLAST  WAVE 


DEFLECTED  BOW  WAVE 


M. 


BOW  WAVE 


EXPANSION  FAN 


NEW  STEADY 
STATE  BOW  SHOCK 


TRANSMITTED  WAVE 


CONTACT  SURFACE 


BASIC  GEOMETRY  AND  NOMENCLATURE  OF  g-D  ANALYSIS 
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FIG. 2. 


PEAK  REFLECTION  PRESSU  RE/PRE-I  NTER  ACTION  SURFACE  PRESSURE,  Frefv  A 


McDONNELL-DOUGLAS  PREDICTION  OF  MACH  REFLECTION 

BOUNDS  FOR  DNA  SLED  TEST. 
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FI  0.3. 


REFLECTION  PRESS URE/4MBI!:NT  PRESSURE 


COMPARISON  OF  CONE  AND  WEDGE  SOLUTIONS 
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FIG.4. 


REFLECTION  PRESSURE/ AM  Bl  ENT  PRESSURE 


COMPARISON  OF  WEDGE  SHOCK-SHOCK  METHODS 


FIG.5. 


DNA  SLED  TEST  4B-B3-COM  PARISON  OF  SEVERAL 

THEORIES  ANDEXPERIMENT. 
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F  IG.6 


PEAK  PRESSURE  BEHIND  THE  REELECTED  SHCCK  WAVE/IN  ITIA.L CONE  PRESSURE 


SMOCK  WAVE  MACH  NUMBER 


FIG.  7  RATIO  OF  THE  PEAK  PRESSURE  BEHIND  THE  REFLECTED  SHOCK 
WAVE  TO  THE  INITIAL  CONE  PRESSURE  VS.THE  SHOCK  WAVE  MACH 
NUMBER- PRIMUS  AND  EXPERIMENT,  9°C0NE  MODELS 
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PEAK  PRESSURE  BEHINOTHE  REFLECTED  SHOCK  WAVE/lNITIAt  CON  E  PRESSURE 


SHOCK  WAVE  MACH  NUMBER 


FIG.  8  RATIO  OF  THF  PEAK  PRESSURE  BEHINOTHE  REFLECTEC  SHOCKWAVE 
TO  THE  INITIAL  CONE  PRESSURE  VS  .THE  SHOCK  WAVE  MACHNUMBE.R- 
.  PRIMUS  AND  EXPERIMENT,  |5°CONES 
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FIG.  9  RATIO  OF  THE  PEAK  PRESSURE  BEHIND  THE  REFLECTED  SHOCK  WAVE  T  O 
THE  INITIAL  CONE  PRESSURE  VS.  MODEL  MACH  NUMBER-  PRIMUS  NOLI 
EXPERIMENT*  9*  CONE  MODELS  * 
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PEAK  PRESSURE  BEHIND  THE  REFLECTED  SHOCK  WAVE/lNITIAL CONE  PRESSURE 


FIGJO.  RATIO  OF  THE  PEAK  PRESSURE  BEMINDTHEREFUECTED  SHOCK  WAVETO 

THE  INITIAL  CONE  PRESSURE  VS. MODEL  MACH  NUMBER-  PRIMUS,,  MOL, 
£  EXPERIMENT,  15°  CONE  MODELS 
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PEAK  PRESSURE  BEHIND  THE  REFLECTED  SHOCK  WAVE/lNITIALCONE  PRESSURE 


SHOCK  WAVE  MACH  NUMBER 


FIG.II  RATIO  OF  THE  PEAK  PRESSURE  BEHIND  THEREFLECTED  SHOCK 
WAVE  TO  THE  INITIALCONE  PRESSURE  VS. THE  SHOCKWAVE  MACH 
NUMBER-THEORYAND  EX  PERI  WENT,  9° CONE  MODELS 
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SHOCK  WAVE  MACH  HUMBER 


FIG.  12  RATIO  OF  THE  PEAK  PRESSURE  BEHIND  THE  REFLECTED  SHOCK  WAVE 
TO  THE  INITIAL  CONE  PRESSURE  VS. THE  SHOCK  WAVE  MACH  NUMBER- 
THEORY  ANDEXPER\MENT,  15*  CONE  MODELS 
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REGULAR  REFLECTION  ON  CONES-PRI MUS  RESULTS 

0c*H.ao,Mb  =  l.a34 


j  I - 1 - 1 - 1 _ I _ l _ i _ t 
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FIG.I3- 


COMPARISON  OF  MAXIMUM  AND  SIDE-ON  PR ESSURE S 

FOR  REGULAR  REFLECTION 


©c«ll^° 

Mb=l^34 


FLIGHT  MACH  NO.,Moa 
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FIG.I4. 


REFLECTION  PRESSURE/AMBIENT  PRESSURE 


COMPARISON  OF  TWO  WEAK  BLAST  CASES  FOR 

VALIDITY  OF  HYPERSONIC  SIMILITUDE 


Mb  =  l.234-  KN  =  Moo  SI  N  9C  =0-9712 

-  Mco  =  5  ,ec=ii.2° 


INTERCEPT  ANGLER, DEO. 


FIG. 15. 
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